Quantum limit for measurement of a weak classical force coupled to a noisy 

quantum-mechanical oscillator 



C. L. Latune, B. M. Escher, R. L. de Matos Filho, and L. DavidovictQ 

Instituto de Fisica, Universidade Federal do Rio de Janeiro, 21941-972 Rio de Janeiro, Brazil 

(Dated: October 12, 2012) 

Precise measurements of tiny forces and displacements play an important role in science and 
technology. The precision of recent experiments, while beginning to reach the limits imposed by 
quantum mechanics, is necessarily spoiled by the unavoidable influence of noise. Here we obtain a 
quantum limit for the uncertainty in the estimation of a resonant classical force acting on a noisy 
quantum-mechanical oscillator. We determine the best initial state of the oscillator and the best 
measurement procedure, thus getting a rigorous and useful benchmark for experiments aiming to 
detect extremely small forces. 
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Several areas of science and technology rely on the ca- 
pacity of measuring tiny forces and displacements [THS] . 
The sensitivity of recent experiments has attained ex- 
tremely high levels [9l412j. so that the limits imposed by 
quantum mechanics start playing an important role. The 
measurements are affected however by the unavoidable 
interaction between these systems and the environment. 
Evaluating this effect is a difficult task, since the deter- 
mination of the ultimate precision limit in the presence of 
noise is still a challenging problem in quantum mechanics. 
Here we determine the quantum limit for the uncertainty 
in the estimation of a weak resonant classical force acting 
on a noisy quantum-mechanical oscillator, which is used 
as a probe for the measurement of the force, through a 
sequence of discrete measurements. Our solution leads 
to the best probe state, for a given average energy of the 
oscillator, and to a measurement procedure that attains, 
in the asymptotic high-energy limit, the so-called "poten- 
tial sensitivity," which defines a level of sensitivity in the 
estimation of a force in the presence of thermal noise that 
cannot be surpassed by any measurement strategy [5] . As 
a specific application, we derive a precise lower bound for 
the uncertainty in the estimation of the force acting on 
a trapped ion under realistic experimental noise condi- 
tions [T^. 

The detection of a weak classical force is an example 
of the general problem of parameter estimation. A typ- 
ical procedure consists in sending a probe in a known 
initial state through some parameter-dependent physi- 
cal process, and measuring the final state of the probe, 
estimating then from this measurement the value of the 
parameter. The precision generally depends on the initial 
state of the probe, on the measurement and estimation 
procedure, on the dynamical process, and on the amount 
of resources used in the measurement (quantified for in- 
stance by the number of probes or the energy of each 
probe). It is given by the Cramer-Rao bound [13 E], 
which relates the uncertainty Sx in the estimation of a 
parameter x to the Fisher information [15] , T{x) , defined 
in terms of the conditional probability density p{^\x) of 



getting the outcome ^ of the measurement when the value 
of the parameter is x: 

Sx^^J ((Xcst - Xtrnef) > j ^, ^ , (1) 

^JvJ^{Xtrne) 



where 



dx 



(2) 



Xcst is the estimated value of the parameter for a possible 
measurement result, xtmc is the true value of the param- 
eter, V is the number of repetitions of the experiment and 
the average in ([T]) is taken over all possible measurement 
results. The above expression holds for unbiased estima- 
tors, for which {xest) = x. In general, the lower bound in 
([!]) is tight for v oo, however if p{£,\x) is Gaussian the 
bound is attainable for any v >\ [T51 - [T5] . Better preci- 
sion is obtained upon increasing the Fisher information 
or the number of repetitions. 

Quantum mechanics imposes restrictions on the preci- 
sion of the estimation, since two outgoing states corre- 
sponding to two different values of the parameter are not 
necessarily distinguishable, and measurements must con- 
form to quantum constraints. On the other hand, quan- 
tum features, like entanglement and squeezing, help to in- 
crease the estimation accuracy beyond the standard lim- 
its |16H23) , yielding better precision for the same amount 
of resources. 

The application of the Cramer-Rao bound to quantum 
theory was initiated by Helstrom [24] and Holevo |25) . 
Braunstein and Caves |26j showed that maximization of 
the Fisher information in (|2]) over all possible measure- 
ment procedures, leading to the so-called quantum Fisher 
information, yields, in the noiseless simple expres- 

sion for the corresponding uncertainty (fl]), which has 
been applied to many different systems Jt] [22l [27] . If 
the initial state of the probe is a pure state \tpQ) and the 
state of the outgoing probe is \ip{x)) = t)'(a;)|V'o)j where 
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U{x) is an a;-clependent unitary operator, then the cor- 
responding quantum Fisher information is four times the 
variance, calculated in the state \ipo), of the operator 



that is. 



n{x) ■.= i j^U{x), 

ax 



J-Q(a;)=4(^o|[A?i(x)]2|V'o) 



(3) 



(4) 



However, the estimation of parameters in the pres- 
ence of noise poses formidable challenges. Only for very 
special situations it is possible to derive analytic lower 
bounds for the uncertainty [28H32] . 

A practical procedure for evaluating the precision of 
estimation was introduced in [331 IM] and further devel- 
oped in [3S]. Here we show that this method leads to 
a tight analytical bound for the uncertainty in the es- 
timation of a resonant classical force acting on a noisy 
quantum oscillator with a fixed average energy, consid- 
ered as a resource. 

We start by considering momentum measurements on 
the oscillator, for an initial Gaussian state. We show next 
that the corresponding Fisher information actually coin- 
cides with the quantum Fisher information for estimation 
of the force, with the same initial states. Furthermore, 
under the constraint of fixed initial average energy of 
the harmonic oscillator, we show that squeezed Gaussian 
states maximize the quantum Fisher information, thus 
yielding the ultimate precision limit for the estimation of 
the force. 

Momentum measurements on the noisy harmonic os- 
cillator. The Hamiltonian of a harmonic oscillator under 
the action of a resonant classical force is given, in terms 
of dimensionless variables, by Hs/hw = {P^ + X^)2 + 
F cos{ujt)X, where, in terms of the momentum p, the po- 
sition X, the mass m, the oscillation frequency w, and the 
force amplitude /, the dimensionless variables are defined 
by P = p/Vmhuj, X = xy/mLj/H, so that [X, P] = i, and 
F — f / \/hjnu^ . The aim here is the estimation of /. We 
note that the resonance condition is the most favorable 
situation for estimating the force. 

Setting X = {a) + a)/^/2 and P = i{a^ - a)/^/2, 
so that [a,a^] = 1, we have, in the interaction pic- 
ture, Hi{t) = hwF cos{Lut){ae~''^* + a1'e'"*)/v^, with 
Ho = hu!/2{P^ + X^) = huj (fltfl +^). In the rotating- 
wave approximation, valid for 3> 1, one neglects oscil- 
lating terms, getting then Hi{t) — hidFX/2, correspond- 
ing to a nunitary evolution implemented by the momen- 
tum displacement operator U{t) = exp{—iujtFX/2). 

Physical insight into the dynamics of the forced har- 
monic oscillator in the presence of thermal noise can 
be obtained from the corresponding Heisenberg-Langevin 
equation for the momentum operator: 



where 7 is the friction coefficient and fjlt) is a Hermitian 
fluctuation force, with (/^(t)) = and {/■y{t)f'y{t')) = 
7(nT + l/2)(5(t — t'), where ut is the average number of 
thermal excitations at temperature T. Integration of ([5| 
yields 

P{t) = P(0)e-T*/2 ^ -p(^) ^ /■* dt7^(i')e^(*'-*)/2 , (6) 



where P(t) = {ujF/j)[1 — exp(— 7t/2)] is the effective 
momentum displacement produced by the force F. This 
shows that, in the presence of friction, the momentum 
displacement due to the applied force is attenuated by 
2[1 — exp(— 7i/2)]/7t with respect to the noiseless dis- 
placement ujFt/2, obtained when 7 — ?> 0. 

From (|6|, the evolution of the momentum variance is 

((AP)2), ^ ry((AF)2)„ + {2nT + 1)(1 - 7y)/2 , (7) 

where the index t stands for the value of the variance at 
time t, and rj = exp(— 7i). 

If the initial state of the harmonic oscillator is Gaus- 
sian, so is the final state pt, since the oscillator is in- 
teracting with a thermal reservoir and the displacement 
is a Gaussian operation. Then the Fisher information 
for the force estimation corresponding to a momentum 
measurement on the final state of the oscillator is 



dp- 



1 



fd{P\pt\P) 



P\pt\P) V 9F 



((AP)2) 



(8) 

where D{ri) — (a;/7)(l — y^) and \P) is an eigenstante 
of the momentum operator. Eqs. ([t]) and ([8| imply that, 
for an initial minimum-uncertainty state in X and P, 



TpiF)^D^{rj) 



4((Al)2>o 



7? + 2(2nT + l)(l-r?)((AX)2)o 



(9) 



dP/dt = a;F/2 - 7P/2 + f^{t) , 



(5) 



where ((AX)2)o = l/[4((AP)2)o] is the variance of X in 
the initial state. 

We show now that this expression coincides with the 
quantum Fisher information for the class of initial states 
considered above, and that no other class of states yields 
a larger quantum Fisher information for the estimation of 
the force. This is accomplished by applying the method 
proposed in Refs. [33H35] . There it was shown that, by 
transforming the non-unitary evolution of a noisy system 
into a unitary one in an extended space, it is possible to 
find an upper bound to the quantum Fisher information, 
which corresponds to the quantum Fisher information 
for the extended system (system plus environment). The 
idea is that measurements on both the system and the 
environment should not lead to less information about 
the parameter than measurements on the system alone. 
It was shown that this bound can actually be attained, 
for some family of extensions. In this case, measurements 
on system plus environment yield redundant information 
with respect to measurements on the system alone. 
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FIG. 1: Beam-splitter model for the coupling with the en- 
vironment. The incoming beam, in mode a, corresponds to 
the harmonic oscillator. Beam splitter Bi induces photon 
losses in the incoming beam, as it deflects photons into mode 
b, which corresponds to the environment. Beam-splitter B2, 
with transmissivity T — > 0, is used to displace the field in 
mode b, through the injection of a high-intensity coherent 
state with amplitude a, such that the product aVT is finite. 
Beam splitter ,83 is used, analogously to B2, to displace the 
field in the environment, upon injection of a coherent state 
with amplitude /3. Beam splitter B4 decouples modes a and 
b, allowing individual measurements on each mode. 



We discuss here in detail the case T ~ 0, which admits 
a simple physical picture. For T 7^ 0, the discussion fol- 
lows similar lines, and is worked out in the Supplemental 
Material. 

Quantum Fisher information for the noisy oscillator. 
Let us consider the harmonic oscillator as a mode of an 
electromagnetic field. The dissipation due to the inter- 
action with the environment R is modeled by the beam- 
splitter Bi, shown in Fig. [T] This device deflects pho- 
tons into mode b (which plays the role of R). If the 
transmissivity of this beam-splitter is 77 = exp(— 7^), then 
the evolution of the variances of the quadratures of the 
electromagnetic field is precisely that given by ([7]), when 
ut = 0, which motivates this beam-splitter picture of 
the dissipation process. The displacement is induced by 
a coherent state sent into a second beam splitter i?2, with 
transmissivity T going to zero at the same time that the 
amplitude a of the coherent state goes to infinity, the 
product a remaining finite. In the Supplemental Ma- 
terial, it is shown that these two operations yield an evo- 
lution for S alone equivalent to the one derived from the 
master equation. The initial state of -I- i? is |V'o)s|0)_r, 
the environment (mode b) being initially in the vacuum 
state. 

After these two operations, the two-mode state 
becomes \i;(F,t))sR = UsR{F,t)\^o)s\0)R = 
e-^FD{n)XsB^\iljo)g\0)j^, where = ^e,{t){aP~a^b} 

is the beam-splitter operator acting on modes a and b, 
with cos 01 (t) = ^Jf\\ a and b are the annihilation opera- 
tors corresponding to modes a and b, respectively; Xs 
is a quadrature operator for mode a, which corresponds 
to the position operator for the harmonic oscillator. 
The Xg-dependent exponential in the above equation 



displaces the momentum as in ([6|. 

In general, measurements on this extended system 
yield more information on the parameter than measure- 
ments on iS* alone. In order to reduce the non-redundant 
information about F in S + R we displace the field in 
mode b, along the same quadrature as the one displaced 
in mode a. This is implemented, as shown in Fig. [T] 
by sending a coherent state with amplitude /3 on beam- 
splitter B^, which has a vanishingly small transmittance, 
as was the case for B2. This is a local operation on R, 
which does not affect S. The evolution operator m S + R 
takes then the form: 



UsR{G,t) 



_ tFGir,}Xn-iFD{v)Xs 



Si 



(10) 



where Xj^ is the position quadrature corresponding to 
mode b. 

Inserting (10) into ^ and Q, we get the respective 
quantum Fisher information: 

J-|«(G) = [-G{v)V^ + Di,j)^]H{{AXsr)o 

+ [Dirj),/T^^ + G{fj)^]H{iAXRf)o,in) 

where the averages are calculated in the initial state of 
S + R. According to the above discussion, this should 
be an upper bound for the quantum Fisher information 
associated to S alone for any value of G. The best up- 
per bound is obtained by determining the function G{ri) 
that minimizes (11). This is done in the Supplemental 
Material. 

The corresponding minimum quantum Fisher informa- 
tion for S + R coincides then precisely with ^ for a zero- 



temperature reservoir, with ((AX)^ 



The generalization of the above procedure for T 7^ 0, 
involving the addition of a squeezing transformation on 
mode a and another environment mode c, yields precisely 
Eq. ([9| - see Supplemental Material. This shows that the 
upper bound obtained by this minimization procedure is 
actually attained for initial minimum-uncertainty states 
in X and P, implying that it is the quantum Fisher infor- 
mation of S alone, for these initial states. It also implies 
that, for these states, the best measurement for the esti- 
mation of the force is a momentum measurement. 

Using ([1]) and the upper bound obtained by minimizing 



( 11 ), one gets, for any initial state of the harmonic oscil- 
lator with variance in X equal to ((AA")'^)o, the following 
inequality for the uncertainty SF in the estimation of the 
force: 



SF > 



1 



'(l-r/)(2nT + l) + 



ri/2 



{{Axy 



(12) 



this limit being attainable, for any integer value of 
v, by Gaussian minimum- uncertainty states in X and 
P. The standard limit corresponds to {{AX)'^)q = 
1/2, so that the right-hand side of Eq. (12) becomes 
[D{ri)V2iy]-'^[2nT{l - v) + ll^^^- 
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FIG. 2: Ratio between the energy- and temperature- 
dependent lower bound for the estimation of the force, ob- 
tained by optimizing the probing time, and the strength 
of the thermal fluctuation, Sf^in, as a function of £ = 
2{{AXf)o{2nT + l). 



Since ( 12 ) is a monotonous function of the variance 



{(AX) )o, one must maximize this quantity in order to 
minimize the lower bound on SF. For fixed average initial 
energy, this is achieved by a momentum-squeezed ground 
state of the harmonic oscillator, as shown in the Supple- 
mental Material, for which {{AX)'^)o ^E+y^E^- 1/4, 
where E is the average energy in units of hu. Since the 
dimensioneless force F is related to the actual force / by 
F = f/Vh muj^, this implies that 



Sf> 




(1 - r]){2nT + I) + 



7?/2 



E 



(13) 

The equality is attained through momentum measure- 
ments on the oscillator, initialized in a momentum- 
squeezed ground state, and using the maximum- 
likelihood estimator [TBHTS] . 

The lossless case is obtained for 7 -> 0; then, 
only the second term inside the brackets is left. 
On the other hand, when E + ^/W^^^TJA > 
r]/\{l — rj){2nT + 1)] the lower bound becomes 



^]~V(l-^)(l + 2nT). Therefore, for 
sufficiently large energy, one gets an expression simi- 
lar to the one for a coherent state (standard limit), no 
matter how small are the losses, but lowered by a fac- 
tor -^7(1 — r]){2nT + l)/[2nT{l — 77) -|- 1], which becomes 
•y/1 — rj when nx <C 1. 

Sequential measurement procedure. If the harmonic 
oscillator senses the force during a time t such that 
7^ ^ 1 and then is submitted to a single measure- 
ment, it follows from Eq. (13) that the uncertainty in 



the estimation of the force will be limited by Sfoo = 
^/m^UJJJ'^cotgh[huJ/(2kBT)]/2, which is independent of 
the total probing time. This is due to the fact that the 
steady-state of the forced probe oscillator, which is in- 
teracting with a thermal reservoir, is independent of its 
initial state. It suggests that, if the unknown force acts 
during a time ttot, it is better to probe the force for 
an appropriate time r, measure the probe system, re- 



set the system and repeat this procedure ;/ times, with 
1^ = ttot/T. With this approach, the total probing time 
is divided into v parts and the uncertainty is limited by 



Sf > 5U 



(1 - l/£)e- 



(14) 



where £ = 2{{AXf)oi2nT + 1). 

The optimal probing time Topt leads to the minimum 
lower bound in Eq. (14 1, denoted by Sfs. Then, for 
£ ^ 1, the strategy of sequential and discrete mea- 
surements, performed in the optimal dimensionless time 
7ropt ^ [£/24]-i/3 ^ leads to 

Sfs > <5/mi„[l + (l/8)(3/£)2/3 + 0{£-')] , (15) 



where (^/min = a/ (2rn/ia;7/itot)cotgh [hio /2kBT) is the 
probe "potential sensitivity" derived in Ref. [3], which 
measures the strength of the thermal fluctuation force 
acting on the oscillator and defines the maximum level of 
sensitivity in the estimation of /, valid for any measure- 
ment strategy. We have thus proven here that this bound 
is actually attainable asymptotically by a minimum- 
uncertainty squeezed state. 

The behavior of Sfs as a function of £ is displayed 
in Fig. |2] for fixed temperature, the gain obtained by 
squeezing the initial state is at most of the order of 10%. 

These results can also be applied to the diffusive situ- 
ation considered in Ref. |1^, by letting 7 —7> 0, tit —J' 00, 
with ^riT = (constant). One gets then from (14) that 



Sf > ,/AmtujV/ttot^l + [4I?((AX)2)ottot]-i • (16) 

In this diffusion limit, the optimal measurement time is 
the total time, which is a consequence of the fact that, 
while the displacement grows linearly with time, the noise 
in the momentum measurement increases with the square 
root ot time. The first factor on the right-hand side 
of the above equation corresponds precisely to the ex- 
pression derived through heuristic arguments in |10| . It 
is seen to overestimate the reachable precision, unless 
4I?((Ai')2)oitot > 1- For the conditions assumed in [10], 
this corresponds to ^tot 3> 1 ms. 

In conclusion, we have been able to completely solve a 
precision problem involving the estimation of the ampli- 
tude of a resonant classical force through measurements 
on a probe consisting of a noisy quantum-mechanical os- 
cillator. The force is estimated through a discrete se- 
quence of measurements on the oscillator, at optimal 
time intervals. We have determined the ultimate pre- 
cision limit, as a function of the average energy of the 
oscillator, and also the best probe state and the best mea- 
surement procedure, thus yielding a rigorous and useful 
benchmark for experiments that aim to detect extremely 
small forces and displacements. 
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SUPPLEMENTAL MATERIAL 
Evolution of the forced harmonic oscillator under a zero-temperature reservoir 

In this Section we derive some useful results concerning the evolution of the forced harmonic oscillator interacting 
with a zero-temperature reservoir. 

Two-step decomposition of the evolution 

Here we show that the simultaneous actions of the external force and the noise can be decomposed into two successive 
operations. The first one being a purely dissipative evolution, corresponding to an interaction of the system with a 
zero-temperature thermal reservoir, and the second one corresponding to a displacement in phase-space. 

In the Markov limit, the master equation that describes the evolution of a harmonic oscillator in the presence of a 
resonant force and a zero-temperature reservoir, in the interaction picture and under rotating-wave approximation, is 
given by 



dp{t) _ .ojF 
dt ~ ~^~2' 



- ^ [a^ dp{t) + pit)^^ a - 2ap{t)a^ . (17) 



The unitary contribution on the right-hand side can be eliminated by choosing a convenient picture, where the density 
operator becomes 

pj,{t) = e^^W^^p(i)e-'^(*)^^ , (18) 
with D{t) being a time-dependent function to be determined. In this picture, the master equation becomes 

= iF - + |£>(i)) [X, Mt)] - I [&^apnit) + PD{t)d^a - 2dpD{t)a^] • (19) 

The function D{t) may be chosen as a solution of the equation 

^ + |i?(t)-a;/2 = 0, (20) 

with the initial condition -D(O) = 0, which leads to the solution 

D{t) = -{l-e-^'/^) , (21) 
7 

noted D{ri) = — (1 — ^) in the following, with 77 = e"'*'*. 

Therefore, in this picture, and for the chosen D{rj), the master equation that describes the evolution of poit) is 

^R^ = -l[a)apD{t)+PD{t)a^a-2dpD{t)a^] , (22) 

which only contains the effects of the interaction between the oscillator and the zero-temperature reservoir. 

This implies that the evolution of the system under the joint influence of the applied force and the interaction with 
the environment can be decomposed into a first step, involving only a purely dissipative evolution, followed by an 
effective 7-dependent phase-space displacement. 

Equivalence between the extended unitary evolution and the master equation treatment 



We show now that there is a convenient purification of the non-unitary evolution stemming from the above master 
equation, which yields a tight bound for the quantum Fisher information corresponding to the estimation of the force 
acting on the harmonic oscillator. 
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The solution of (22) may be found, for instance, in Ref. 



Going back to the interaction picture, one can see 



that the harmonic oscillator will evolve as 



oo 



^tD{-n)FX 



(23) 



where po = \ipo) s s{^o\, and \'4'o)s is the initial state of system S. 

This non-unitary evolution can be seen as resulting from an unitary evolution on an enlarged Hilbert space com- 
prising the system and an environment, when this environment is not monitored. Here, we represent the environment 
by another harmonic oscillator and show that the unitary interaction between the system and this environment is 



enough to implement the non-unitary evolution given in Eq. ( 23 1 
Let then 



n=0 



-iDiri)FXs 



(1 ^)"^^^ata/2~n/ft\n 



(24) 



be a purification of p{t), that is, if one traces out the environment in |^'(t))(^'(i)|, one is left with p{t). In the above 
equation, \n)fi are Fock-states of the environment R, and a (b) is an annihilation operator corresponding to system 
S (environment R). The subindex S has been added to the position operator corresponding to system S. Now, it is 
straightforward to show that \'^{t)) can be rewritten as 



(25) 



where Bi = e^recos (v^7)(afc^-atb) 

can be seen as the transformation performed by a beam-splitter with transmissivity 
rj on the input modes represented by a and b. Since |5'(i)) is a purification of p{t), the evolution described by the 



master equation (17) can equivalently be modeled by a beam-slpitter-like unitary interaction between the system S 



and the environment R (represented here by a harmonic oscillator) , followed by a displacement in phase-space of the 
system S, when the environment R is not monitored (is traced out). 



Minimization of the quantum Fisher information of system plus environment 
The quantum Fisher information corresponding to system S plus environment R is given by 



J-|^(G) - [-G{rj)^/T^ + D{fj)^]H{iAXsr)o + [D{r,)^T^+G{v)^]H{{AXn)\. (26) 
Minimization of this quantity over G leads to 



(1 - v){{^Xsf)o + v{{^XRf)o\ G{v) + [{{^Xnf)o - {{AXs)\\ D{rj) = 0. 

The minimum is, therefore, reached for 

yrKl^ [((AX5)2>o - ((Al^)2)o" 



Gopt(?7) = 



D{rj). 



(27) 



Inserting this quantity into the expression of Tq^{G) and using ((AXfl)^)o = 5, since the environment R is initially 
in the vacuum state, we obtain 



J-|«(Gopt)-i?^(^) 



4((AXs)2)o 



r; + 2(l-,?)((AXs)2)o 



Evolution of the forced harmonic oscillator under a thermal reservoir with T 7^ 



In this Section, we generalize the previous results to the case of non-zero temperatures. 
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Two-step decomposition of the evolution 

Here we show that the simuhaneous actions of the external force and the thermal noise can also be decomposed 
into two successive operations. 

When the probing harmonic oscillator interacts with a thermal reservoir at a temperature T, its evolution, in the 
interaction picture, is dictated by the master equation 



dm 

dt 



-i^ X,p{t) - I (riT + 1) [a^a/3(t) + /5(t)a^a - 2ap(t)a'f] 
[aa}p{t) + p(t)aa'f - 2a} p(t)a\ , 



(28) 



where is the mean number of thermal excitations in the reservoir, which obeys the Bose-Einstcin distribution, 
given by 



1 



riT 



ahu)/kBT 2 



(29) 



whith fcs being the Boltzmann constant. 

This master equation is simplified by using the same procedure as in the zero-temperature case, that is, one defines 

pc(t) = e'^('')^^p(t)e-»^('')^^ , (30) 

with 13(77) defined as before. One finds then that pD satisfies the above master equation without the force term. 

Equivalence between the extended unitary evolution and the Master equation treatment 



The solution of ( 28 ) without the force term may be found in Ref. 



(riT + l)'(?^T)■'^'l ,.t\5-i- r-r\i-i 
l\j\r2 



(31) 



where the functions ri, r2 and are defined by 



l-rj 



K(l - r?) + 1] 

nril — 77) + 1 
1 



717^(1 — 77) + 1 

and 77 = e""*"*. Therefore, in the interaction picture, the evolution of p{t) is 

1+3 



p(i)=e-^^('')^^<(r3ei"['-=l'^''^ ^ 



l\]\r2 



(32) 



A purification of p{t) can be built from this solution. Notice that, to find this purification, one considers an 
environment with two harmonic oscillators. Such a purification may be given by 



,-.z,(.).x.^>nMata g J (nr + l)'(y)^r^V3 (d^i^V ^^\^,)s\OUAO)n. 



1,3=0 



j! l\ 



where the states \l) and |j)i?2 are Fock states of the environments Ri and R2 respectively, and 6 (c) is the anni- 
hilation operator for the environment i?i (i?2)- The subindex S has been added to denote an operator acting on S. 
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This purification may be seen in terms of three unitary evolutions: the first one corresponds to a beam-sphtter-hke 
interaction between the system S and the environment the second one corresponds to a two-mode squeezing- hke 
interaction between the system S and the environment i?2 , and the third one corresponds to phase-space displacement 
in S space. Defining the first two unitary operators by 

S = ge2(t)(atct-ac) 

where Oi(t) and 92(t) are given by 



V 



nril - 1]) + 1 



9i(t) — arccos 
h{t) = arccosh \/nT{l — rj) + 1 , 



then, the above purification can be rewritten as 



(33) 



Minimization of the quantum Fisher information of system plus environment 



One should note that there are many purifications of the state p{t) of S, corresponding to different states of 
S + Ri + i?2- Indeed, for any Hermitean operator Hi 2 acting only on the environments -Ri and R2, the pure state 
1$), given by 



1$) = e-*-"('')^^i-^|«'(i)) , 



(34) 



is another purification of p{x). Notice also that Hi 2 should be chosen properly in order to minimize the quantum 
Fisher information in S* + Ri + R2 ■ 

An upper bound to the quantum Fisher information in S may be calculated with 1$). This bound is given by 



l""-^- = [2i?(,7)]\(0|^^(0| <,(^o| {a \b\S^ (Is + ^1,2) ^i^i] }' ^^)s\G)rMr. 



(35) 



A possible choice of the operator Hi 2 aimed to erase part of the non-redundant information about the value of the 
force F in |$), as compared with the same information in /5(t), is Hi;2 = AiX/j^ -|- \2X^^^ where X/j^ (^^2) 
the position operator of the oscillator in Ri {R2) space. This choice is based on physical insights on the enlarged 
unitary process. Since |$) is an entangled state in S + Ri + R2, the effect of a phase-space displacement along the 
momentum-axis in S may affect the state in 5* -I- i?i -I- i?2 , by changing the correlations between the system and the 
environments. Indeed, for Hi 2 — 0, after disentangling S + Ri + R2 with the operation bIS^ , which does not change 
the quantum Fisher information, the effective unitary evolution in S + Ri + R2 is 



s,Ri,R2 = exp' 



\^-iD{7])F cosh [02(0] cos [Oiit)] Xs - cosh [02(0] sin [6*1(0] ^Ri + sinh [02(0] ^R2 } 



(36) 



It is clear that, with Hi 2 = XiX^-^ + A2A^i?2 ^^d convenient values of Ai and A2, it is possible to erase at least part 

of the non-redundant information in |^'(t)).^ 

After a straightforward calculation, Eq. (351, with the above choice of Hi^2, can be rewritten as 

J-^«^«^(Ai,A2) 



- ((AX5)^)o 

+ ((AIhJ^)o 
+ ((AXhJ2)^ 



(cosh[02(O] + ^2 sinh[02(O]) cos[0i(O] + -^i sin[0i(O] 



1 2 



-(cosh[02(O] + A2 sinh[02(O]) sin[0i(O] + Ai cos[0i(O] 



sinh[02(O] + A2 cosh[02(O] 



(37) 
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Then, the optimal vahies of Ai and A2, which minimize the above equation, are 



(opt) _ 



(((Alflj2)o - ((AXs)2)o)cos[gi(t)]sin[gi(t)] 
A^°P*' = -cosh[02(O]sinh[02(t)] 

{{AXs)')o{{AXnJ')o 



(cosh[02(t)] + A^°P*)sinh[^2(t)]) 



{{AXnJ')o + 



cosh^ [02 (t)] {{AXR^r)o + sinh^ [e^ (t)] 



{{AXB.,rhcos^Oi{t)] + {{AXsrhsin^Mt)] 
{iAXsr)o{{AXnJ^)o 



{{AXn,y)ocos^[Oi{t)] + {{AXsy)osm'Mt)] 



Therefore, the minimum value of J-^^i^^(Ai,A2) is given by 



sinh^ [02 (t)] cosh^ [02 (t)] sin^ [0^ (t)] cosh^ [02 (t)] cos^ [0^ (t)] 



((AXhJ^>o 



{{AX- ^2 



B.2) /o 



{{AXr^^ 



(38) 



Substituting 0i{t) and 02{t) in terms of rj and nx, and setting {{AX]^-^)'^)q = {{AX]i^)'^)q = 1/2 in the above equation, 
since both environments start in a vacuum state of the harmonic oscillator, we get the following upper bound 



_^Sflifl.(;^(opt)^^(opt)^ = [2D{7^)] 



-\ -1 



2(l-r;)(2nT + l) + 



(39) 



Maximization of {{AXf)o 

We want to maximize the variance {{AX)'^)q under the constraint of fixed resource {X'^)q + (P^)o = and taking 
into account the physical restriction imposed by the Hciscnberg uncertainty relation {{AX)'^)q{{AP)'^)q = a > 1/4. 
This leads to a system of two equations with two unknown parameters: 

{{AX)\{{AP)\ = a (40) 
{{AX)\ + {{APf)o = 2E* (41) 

with 2E* ^2E- {X)l - (P)o- 

The solutions are ((AX)2)± ^ E* ± V£^*2 - a. It follows that the maximum value of the variance in X under the 
two constraints imposed above is ((AX)2)q = E* + ^/E*"^ — 1/4, reached when a — 1/4, that is, when the state is a 
minimum uncertainty state in X and P. 

Minimum-uncertainty states in X and P are pure Gaussian states, which can be obtained from the ground state 
through squeezing characterized by a real parameter r and a complex displacement z. Because the variances ((AAl )2)o 
and ((AP)2)g do not depend on the value of z, the best strategy to maximize the variance {{AX)'^)q, for a fixed value E 
of the average energy, is to invest all the energy in the squeezing of the state. This allows one to conclude that the state 
that maximizes the variance of X, under the constraints imposed above, is a squeezed ground state j-^o) = »§'(?') |0), 
with S{r) = exp §(a2 - a^^) and r = 1/2 In [2{E + y'E^ ~ 1/4)]. 
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